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Abstract. A procedure is described for denning a generalized solution for stochas- 
tic differential equations using the Cameron-Martin version of the Wiener Chaos 
expansion. Existence and uniqueness of this Wiener Chaos solution is established 
for parabolic stochastic PDEs such that both the drift and the diffusion operators 
are of the second order. 



1. Introduction 
Consider a stochastic evolution equation 

(1.1) du{t) = (Au(t) + f(t))dt + (Mu(t) + g(t))dW{t), 

where A and M. are differential operators, and W is a Wiener process on a proba- 
bility space (Q, J 7 , P). Traditionally, this equation was studied under the following 
assumptions: 

• Al. The operator A is elliptic, the order of the operator M. is less than the 
order of A, and A — \M.M* is elliptic (possibly degenerate) operator, 

In fact, it is well known that unless assumption Al holds, equation (11.11) has no 
solutions in L2(fi; X) for any reasonable choice of the state space X. 

It was shown recently (see jH [5], [6] and the references therein) that if only the oper- 
ator A is elliptic and the order of M. is smaller than the order of A, then there 
exists a unique generalized solution of equation (II. ip . This solution is often re- 
ferred to as Wiener Chaos solution. It is given by the Wiener chaos expansion 
u {t) = J2\ a \<oo u a (t) where {^ Q }| Q ,| <00 is the Cameron-Martin orthonormal basis 
in the space L 2 (fl] J rW ; X) of square integrable random elements in X measurable 
with respect to to the sigma-algebra J rW generated by the Wiener process. The 
Cameron-Martin basis is indexed by multiindices a = (ax, ■■■) ■ It was shown 
that for certain positive weights Q = {q (a)}| a i <00 , the weighted norm 



1ii] \q,x : = q2 ( a ) IKHL((0,T);X) < °°> 

|o| <oo 
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where X is the appropriate Hilbert space characterizing the "regularity" of the solu- 
tion. Note that without assumption Al 

E \\ U \\l 2 ((0,T);X) = W Ua (^Wl 2 ((0,T);X) = °°- 

|«|<oo 

In this paper, we consider the Cauchy problem for the following stochastic partial 

differential equation: 

(1.2) 

du = (aijDiDjU + biDiU + cu + f)dt+(pijD i DjU+<TiDiU + b , u+g)dW, t e [0, T], iGl. 

In contrast to the previous work, this is a parabolic SPDE of the full second order, 
in that the drift and diffusion operators have the same order 2. We construct a scale 
of weighted Wiener chaos spaces (related but not identical to Kondratiev's spaces) 
and prove the existence and uniqueness of the solution in the spaces from this scale. 



2. Constructing a Solution: an Example 

Let F = (f2, JF, {jF t } <t<T, P) be a stochastic basis with the usual assumptions and 
W = W(t), < t < T, a standard Wiener process on F. For a Hilbert space X, denote 
by L 2 (W;X) the collection of X- valued random elements that are square integrable 
(E|| • ||^ < oo) and are measurable with respect to the sigma-algebra generated by 

W{t),te[o,T}. 

Consider the Ito equation 

(2.1) u(t,x) = e^ 2/2 + f u xx {s,x)ds+ [ u xx (s, x)dW(s), t 6 [0, T], x e R. 

Jo Jo 



If there is a solution, its Fourier transform in space, u(t, y) = (1/ v27r) f R e txy u(t, x)dx 
satisfies 



(2.2) u(t,y) = e 



-v 2 /2 _ ,.2 



y 2 [ u{s,y)ds-y 2 [ u{s,y)dW{t), te [0,T], ye 
Jo Jo 



For each fixed y, ( 12.2ft defines a geometric Brownian motion: 

(2.3) u(t,y)=e' (1+t ^- ( y 4 ^ t -y 2w(t l 
Let i7 7 (R) be the Sobolev space 

(2.4) (/: / (l + \y\ 2 y\f(y)\ 2 dy<oo\. 



Since 

(2.5) E\u(t,y)\ 2 = e- 2 ^ y2+ y 4 \ 

the solution of <^B) cannot be an element of L 2 (W; L 2 ((0, T);H^{R))) for any 7 e R, 
even though the initial condition is non-random and is an element of iY 7 (IR) for every 
7 G R. 
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Let us try another approach. Once again, assuming that the solution exists, we apply 
the Ito formula to the product u(t,x)£h(t), where 

(2.6) £ h {t) = exp Qf h(s)dW(s) - ^ jf h 2 (s)ds^j 
and h = h(t) is a smooth deterministic function. Since 

(2.7) £ h {t) = l+[ £ h (s)h(s)dW(s), 

Jo 

we conclude that the function 

(2.8) u h (t,x)=E(u(t,x)£ h (t)), 
if defined, must satisfy the heat equation 

(2.9) u h (t, x) = e~* 2 / 2 + j\l + Hs)) d2u ^ x) ds. 

If sup t \ h(t)\ < 1, then this equation has a unique solution in every if 7 (R) and 

(2.10) u h (t, x) = Eexp(-(X(t, x)) 2 /2), 
where 

(2.11) X(t,x)=x+ f ^/2{1 + h{s))dW{s). 

Jo 

In other words, while existence of a solution of equation (12.1 j) is still unclear, we now 
have a family of functions Uh(t, x) defined by (12.101) . All we need now is a systematic 
procedure of relating the family of deterministic functions Uh = Uh(t,x) to a random 
process u = u(t,x); then this process is natural to call a solution of (12. 1|) . 

Here is a possible way of constructing a stochastic process from Uh- Let m = {m^, k > 
1} be the Fourier cosine basis in L 2 ((0, T)): 

(2.12) mi ( s ) = _L m k {t) = ^cos( * {k ~ l)t \k>l- 0<t<T. 



Then 

(2.13) h(t)=^2h k m k (t), 



k>l 



For every fixed t e [0, T] and 7 £ K, we can now interpret the function u^(t, •) as 
a mapping from the set of sequences h = (hi, h 2 , ■ ■ .) to the space iT ) '(R a! ), and, 
as equalities (12.101) and (12.111) suggest, this mapping is analytic in the region {h : 
J2k>i hi < £ } ^ or sufficiently small e. We will now compute the derivatives of this 
mapping. 

Let J be the collection of multi-indices a = {a^, k > 1}. Each a G J has non- 
negative integer elements a k and 

(2.14) \a\ = ^a k < 00. 
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We also use the notation 

(2.15) a! = J]« fc ! 

k 

and consider special multi-indices, a = (0) with \a\ = and a = with |a| = 1, 
ati = 1. 

For each a G J define 

1 d^u h {t,x) 



(2.16) u a {t,x) 



y/ddh^dh?-- 

Then 

(2.17) u h (t,x) = y^ y u a (t,x) — 



h=0 



V Oil 



where 

(2.18) h a = J[h?. 

k>l 

On the other hand, by direct computation, 

(2.19) £ h {t) = E(S h (T)\^) = E ^nUt), 
where 

(2.20) Ut) =E(^|^T), £ a =^=T[H ah ( f m k (t)dW(t)) 



and 



(2.21) if„(x) = (_i)^ 2 / 2 « e -V 



2 



is n-th Hermite polynomial. It is a standard fact PQ that the collection {£ Q , a G J7"} 
is an orthonormal basis in /^(W; R). 

The functions w a (t, x), a G J7", uniquely determine Uh(t,x) according to (12.171) . On 
the other hand, if 

(2-22) ElK(0llW) <oo > 

then the if 7 (IR)-valued random process 

(2.23) u(t,x) = J2 u *(t,x)ta 

satisfies E(u(t,x)£h(t)) = Uh(t,x); if, in addition, u is jF^-adapted, then also 

(2.24) u(t,x) = ^u a (t,x)£ a (t). 

aej 

If condition (12.221) fails, then (12.231) is a formal series, which we define to be the 
stochastic process corresponding to the family Uh- 
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As (12.51) suggests, if Uh is the solution of fl2.Qj) . then (12.221) fails for every 7. Let us 
now see how fast the series diverges. Equality (I2.9P implies 

u (0) (t,x)=e- x /2 + J M = °; 

f* d 2 u e (s,x) , /"* d 2 U( )(s,x) , . , 

(2.25) = y o g; 2 ; w jf ( °;V w *)^, i«i = i ; 

f t d 2 u a (s,x) ^ f 
Jo ox k=i J 



l t d 2 u a (s,x) — f f d 2 u a _ ek (s,x) 
Ua{t) = / ^ ds+y y/Ok — m k (s)ds, \a\ > 1. 



cte 2 

Equations of the type (12.251) have been studied [U Section 6 and References]. In 
particular, it is known that 

(2-26) £ |K(0IIW) = ^P^o||^ (R) , 

|a|=n 

where D x = d/dx, $t is the heat semigroup, and Uo(x) = e~ x2 ^ 2 . To simplify further 
computation, let us assume that 7 = 0. Then, switching to the Fourier transform, 



(2.27) \\Dl^M\l 2m = I \yre-y 2 ^dy 

JR 



r (2w + \) 
(i + t) 2n . 

Using Stirling's formula for the Gamma function T, 

( ? /f\ 2n 

(2.28) IMOHLcr) = (r+t) c ^ n[ > 

\a\=n \ "r / 

where the numbers C (n) are uniformly bounded from above and below. Similar result 
holds in every if 7 (R). Thus, (I2.22p does not hold, but instead, by (12.281) . we have 



(2-29) E AlM*)llW) <0 °- 



We denote by (£) 0) o(W; # 7 (K)) the collection of formal series fl2T24l) satisfying (jgj29j) ; 
the reason for using (£)o.o in the notation will become clear later. Note that we had 
equalities in all computations for equation ( 12. II) that lead to ( 12.291) . which suggests 
that (£)o,o(W; if 7 (IR)) is the natural solution space for equation (12. ip . For a more 
general stochastic parabolic equation of full second order in M. d , the natural solution 
space turns out to be (£) M (W; L 2 ((0, T); # 7 (M d ))) for suitable p, q < 0. 

In the next section we address the following questions: 

(1) How to define the spaces (£) Pi9 (W;X) for p, q 6 R without relying on an 
ortho normal basis in L 2 ((0, T))? 

(2) How to construct a solution of a general stochastic parabolic equations of full 
second order? 
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3. General Constructions and the Main Result 



As before, let F = (ft, J 7 , {J-t}o<t<T, P) be a stochastic basis with the usual assump- 
tions and W = W(t), < t < T, a standard Wiener process on F. Denote by 
H s = H s ((0,T)), s > 0, the Sobolev spaces on (0, T) with norm 

(3- 1 ) II ' lis = ll A * /2 ' \\l 2 ((o,t)), 

where A is the operator 

T 2 d 2 

with Neumann boundary conditions. This norm extends to functions of several vari- 
ables via the tensor product of the spaces H s . 

Definition 3.1. Given real numbers p,q and a Hilbert space X, (£) Pi9 (W;X) is the 

closure of the set of X -valued random elements 

(3.3) 

V = Vo + f2[ [ ■■■[ Vk(si,...,s k )dW(s 1 )...dW(s k - 1 )dW(s k ), N>1, 
k=1 Jo Jo Jo 

with respect to the norm 

(3-4) \\v\\ 2 p>q ;x = \\Vo\\ 2 x + Yl Jkiy H H^H* H*> 

where each rjk, k > 1, is a smooth symmetric function from [0,T] k to X. 
Remark 3.2. (a) It is known [21 Ej that, for rj of the type (13.31) . 

N 



(3-5) nv\\x = \\Vo\\x + J2y\\\\Vk\\o\\x 



fc=i 



(b) The definition of each individual (£) p j3 (W; X) inevitably involves arbitrary 
choices, such as the norm in H q ((0,T)). Further analysis shows that different choices 
result in shifts of the indices p,q, and the space U Pj(? (£) Pig (W; X) does not depend on 
any arbitrary choices. In the white noise setting, where Q is the space S'(Mr) of the 
Schwartz distributions and P is the normalized Gaussian measure on S, the inductive 
limit U Pi9 (£) Pi9 (W; R) is the Kondratiev space (<S)-i [2]- 

Remark 3.3. If X is the Sobolev space i? 7 (M d ), then we denote the norm || • || p , 9i x 

|| ■ llp,<j;7- 

(3-6) || ■ ||p,g;_H'7(Rd) = || • || p ,g; T 

Proposition 3.4. Let rj = f£ a , where f G X and £ a zs defined by (12.201) . T/ien 

(3-7) ll^llLx = ^N 2 ^||/|| 2 ,, 

where 

(3.8) N 2?a = JJP 9Qfe . 

fc>i 
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Proof. Let \a\ — n. It is known [3] that 

(3.9) £, = __/ / . . . / E a {s x , . . . , s n )dW{s x ) . . . dW(s n ^)dW(s n ), 

va! Jo Jo Jo 

where E a is the symmetric function 

(3.10) E a (si, . . . , s n ) = m h (s a ^)) . . .m in {s a {n))- 

<?ev n 

In (I3.10p . the summation is over all permutations of {1, . . . , n}, the functions m k are 
defined in (12.121) . and the positive integer numbers i\ < ii < . . . i n are such that, for 
every sequence (b k , k > 1) of positive numbers, 

(3.H) \{b a k k =b n -b l2 -...-b ln . 

k>l 

For example, if a = (1,0,2,0,0,4,0,0,...), then \a\ = 7 and i\ = 1, i 2 = ^3 = 3, 
14 = . . . = if = 6. Thus, in the notations of (13.41) . we have 

L -.E a f, if k = n, 



(3.12) % =<Va! 

v 0, otherwise 

Note that 

(3.13) \\E a \\ = \fa\Vn\. 

By definition ( 13. 2D of the operator A we have 

(3.14) A q/2 m k = k q m k or \\E a \\ 2 q = N 2qa a\n\. 

The result now follows. □ 



Corollary 3.5. A formal series 
(3.15) »7 = l>a£ 



wrat/i n Q e X ; zs an element of (£) Pj(? (W; X) z/ and on/?/ z/ 

(3-16) > — -- — WVaWx < °°- 

z — ' a ! 



Proof. This follows from (13. 14|) and the equality 

(3.17) \\E a + E p \\l = \\E a \\l + \\Ep\\l a^(3. 



□ 



Denote by (£)™(W) the Hilbert space dual of (£) M (W;R) relative to the inner 
product in L 2 (W;IR), and by ((•,•)) the corresponding duality. In the white noise 
setting, n Pig (£) p ' 9 (W) is the space (S)\ of the Kondratiev test functions [2j. If n G 
(£) M (W; X) and C e (£) M (W), then ((n, C» is defined and belongs to X. 
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For h G L 2 ((0, T)), define 

(3.18) £ h = £ h (T) = exp fj\(s)dW{s) - ±£ \h{s)\ 2 ds\ . 

Proposition 3.6. The random variable £h is an element of (£) p,q (W) if an only if 

(3.19) \\h\\^ q <2 p . 

Proof. Since 

(3.20) £ h (T) = 1 + / h{t)S h {t)dt, 

Jo 

it follows that 

(3.21) £ h = 1 + V / / . . . / h(s k ) ■ ■ ■ h( Sl )dW( Sl ) ■ ■ ■ dW(s k ^)dW(s k ). 

k=1 Jo Jo Jo 

By (E3D and fl3~3|l . £ h G (£) M (W) if and only if 

oo , 

(3.22) J2 { 2 ~ P \\ h W- q ) <°°' 

k=i 

that is, \\h\\ 2 _ q < 2 p . □ 

Definition 3.7. We say that the function h is sufficiently small if (13.191) holds for 
sufficiently large (positive) —p,—q. 

Proposition 3.8. If u E |J q (£) Ptq (W; X) and h is sufficiently small, then 

(3.23) u h =((u,£ h )) 
is an X -valued analytic function of h. 

Proof. For every u G IJ g (£) p , g (W; X), there exist p,q such that u G (£) Pj3 (W;X); 
by Proposition 13. 6\ Uh will indeed be defined for sufficiently small h. Similar to (12.171) 

we have 

(3.24) = E ^TT 

and this power series in /i a converges in some (infinite-dimensional) neighborhood of 
zero. □ 

From now on, Di = d/dxi, and the summation convention is in force: Cjdj = Yli c A^ 
etc. 

Consider the linear equation in M d 

(3.25) du = (a^DiDjU + biDiU + cu + f)dt + (pijDiDjU + o~iDiU + vu + g)dW 
with initial condition u(0,x) = v(x), under the following assumptions: 

BO All coefficients are non-random. 
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Bl The functions a^- = a,ij(t,x), pij = pij(t,x) are measurable and bounded in 
(t,x) by a positive number C , and 

(i) 

\aij(t,x) - aij(t,y)\ + \pij(t,x) - Pij(t,y)\ < C Q \x -y\, x,y eR d , < t < T; 

(ii) the matrix (a^) is uniformly positive definite, that is, there exists a 5 > 
so that, for all vectors y G M. d and all (t, x), a^yiyj > S\y\ 2 . 
B2 The functions bi = bi(t,x), c = c(t,x), o~i = <Ji(t, x), and v = v(t,x) are 
measurable and bounded in (t, x) by the number C . 

B2 

(3.26) u G U(£) Pi(/ (W;L 2 (R d )), f,g G |J(£) M (W; L 2 ((0, T); i/- 1 ^))). 

For simplicity, we introduce the following notations for the differential operators in 

(3.27) A = CLijDiDj + biDi + c, B = p ii D i D j + o^A + v- 

Definition 3.9. A solutionu of f[3~25|) is an element o/U M (£) M (W; L 2 ((0,T); f/" 1 ^))) 
snc/i that, for all sufficiently small h and all t G [0, T], the equality 

(3.28) u h (t,x) =v h (x) + [ (A + h(s)B)u h (s,x)ds 

Jo 

holds in H- x {R d ). 

The following theorem is the main result of this paper. 

Theorem 3.10. Assume that, for some p > and q > 1, uo G (£) Pi <j(W; L 2 (M d )) and 
f,g are elements of the space (£) Pi9 (W; L 2 ((0,T); if _1 (IR d ))) . Then there exist r, £ < 
siic/i t/iat equation (I3.25P nas a unique solutionu G (£) r ^(W; L 2 ((0, T); if 1 (M d ))) and 

(3-29) jT \\u{t)t Al dt < C ■ (\\v\\l q;0 + £ (\\f{t)\\lr,-i + W)llL-i) * 
The number C > depends only on 5, Cq,p, q, r, £, and T. 

Proof. The proof consists of two steps: first, we prove the result for deterministic 
functions v, f, g and then use linearity to extend the result to the general case. 

Step 1. Assume that the functions v G L 2 (IR d ), f,g G L 2 ((0, T); H~ 1 (R d )) are deter- 
ministic. Then Vh = v, fh = f, gh = 9, and classical theory of parabolic equations 
shows that, for sufficiently small h, equation (I3.28P has a unique solution Uh and the 
dependence of Uh on h is analytic. 

As in the previous section, we write 



(3.30) 



u(t,x) = ^2u a (t,x)£ a 
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where the coefficients u a satisfy 

U(o)(t,x) = v(x) + / (Au(p)(s,x) + f(s,x))ds, 



(3.31) u £k (t,x)= / Au €k (s,x)ds + / (Bu(p)(s, x) + g(s, x))m k (s)ds 



/ Au a (s, x)ds + > Jotk / Bu a _, , , 
7o V ^ 



m q (s,x)= / ^4w a (s, s)c?s + > / Bu a ^ ek (s, x)m k (s)ds, \a\ > 1. 



Denote by $ = £ > s > the semigroup generated by the operator A. It follows 
by induction on \a\ that 

(3.32) 

U( )(t,x) = ®t,ov(x) + / $ tjS /(s)ds, 



w £fc (t,x)= / $ t)S (BM( )(s,x) +^(s,x))m fe (s)ds, 
Jo 

u a {t,x) = —=j I ... I $ t>Sn B® SnjSn _ 1 ---B$ S2)Sl (Bu(o)(si,x)+g(si,x)) 
V a' Jo Jo Jo 

E a (si, . . . , s n )dsi . . . ds n , \a\ = n > 1. 

Therefore, using the usual parabolic estimates, 
(3.33) 

IPaWllifl(Rd)^ — a , Mrlli 2 (Mrf) + y (JI/(*)llff"-l(]R d ) + ||Pw|lH- 1 (K d ) J 

and then (|Ogj) follows from fl3~TB|) . 



Irft 



S'tep As in Step 1, existence and uniqueness of solution follows from unique solv- 
ability of the parabolic equation (I3.28p . and it remains to establish (13.291) . 

Denote by u(t,x; V,F,G,j), 7 G J, the solution of (13.251) with v = V£ 7 , / = F£ 7 , 
g = G£ 7 . Ev = Y.a&j v ^a, etc., then 

(3.34) u(t, x) = u (t, x \ v i, fj, 9~i, l)- 

■yej 

It follows from (I3.3ip that u a (t, x; V, F, G, 7) = if |a| < I7I and 

(3 35) W<, x; V, F, G, 7) _ u « (*> ^ v%' v%> v%' ^ 



y(cr+ 7)! 
Using the results of Step 1, 



(3.36) 



T 

\\u(t,-;v 7 ,f^,g 7 )\\r A1 dt 







Now (13.291) follows from (I3.34p by the triangle inequality. □ 
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